Shape coexistence and configuration mixing in the neutron-deficient Hg isotopes 
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The coexistence of prolate and oblate geometrical shapes in the neutron-deficient mercury isotopes 
is analyzed on the basis of the Gogny energy density functional. Level energies and transition 
rates are calculated within the interacting boson model (IBM) with configuration mixing, whose 
parameters are determined by mapping the potential energy surface, obtained from the constrained 
self-consistent mean-field calculation employing the Gogny DIM functional, onto the IBM analog. 
The triaxial degree of freedom turns out to play an important role in the description of the shape 
coexistence and the configuration mixing in the considered mercury nuclei. An optimal choice of 
the configuration mixing IBM Hamiltonian is suggested, that meets the indication of the Gogny- 
D1M energy surface. The collective low-lying structure in the considered mercury chain is described 
nicely without adjustment, in particular the oblate- prolate band structures of the near mid-shell 
nuclei « 4 -i88 Hg . 



PACS numbers: 21.10.Re,21.60.Ev,21.60.Fw,21.60. Jz 



I. INTRODUCTION 

Both shape transition and shape coexistence in finite 
nuclei have been a theme of major interests in the field 
of low-energy nuclear structure study. In some nuclei, 
significantly low-lying excited + states close in energy 
to the + ground-state have been observed, that reveal 
the coexistence of different intrinsic shapes. Numerous 
efforts have already been made to better understand the 
nature of such stunning shape phenomena from both the 
theoretical and experimental sides (see Refs. [U-Ifl for re- 
view) . 

In terms of the nuclear shell model the emer- 

gence of low-lying excited + states can be traced back 
to multiparticlc-multihole excitations. This scenario ap- 
plies to the neutron-deficient nuclei in the Lead region 
with 2 w 82. In this case, two or four protons are excited 
from the Z = 50—82 major shell across the Z = 82 closed 
shell to the /i g / 2 orbit. The residual interaction between 
the valence protons and neutrons becomes subsequently 
enhanced, leading to the lowering of the excited + en- 
ergies. This effect is most significant around the neutron 
mid-shell N = 104. 

In the neutron-deficient Hg isotopes, the weakly de- 
formed oblate ground-state band coexists with a more 
deformed prolate one. Empirically, in these isotopes the 
oblate ground-state band stays constant in energy as a 
function of the neutron number down to 188 Hg, where it 
is crossed by intruder prolate configurations. The band- 
head of the intruder prolate band, i.e., the second + 
state, reaches a minimum in energy around 182 Hg and 
then goes up from 180 Hg [l2|. From the experimental 
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point of view, this picture is quite vividly observed in 
the parabolic trend of the states belonging to the intruder 
prolate band as functions of the neutron number in the 
Hg chain. Given the recent advances in the experimental 
studies on the neutron-deficient Hg isotopes jL3 - [l5l | , it is 
quite timely, as well as significant, to pursue a theoreti- 
cal description of the relevant spectroscopy that can be 
compared with available data. 

Let us stress, that these heavy mass systems are cur- 
rently beyond the reach of large-scale shell-model studies 
with a realistic nucleon-nucleon interaction. Therefore a 
drastic truncation scheme is required to make the prob- 
lem more feasible. Such a framework is provided by the 
interacting boson model (IBM) [l6j], which employs the 
monopole s and quadrupole d bosons, associated to the 
J = + and 2 + collective pairs of valence nucleons, re- 
spectively [l7rll9| . The collective levels and the transi- 
tion rates are generated by the diagonalization of the bo- 
son Hamiltonian composed of only a few essential inter- 
action terms. To handle the configuration mixing in the 
IBM framework, Duval and Barrett proposed to extend 
the boson Hilbert space to the direct sum of the config- 
uration sub-spaces corresponding to the 2np-2nh excita- 
tion (n = 0, 1, 2 . . .) that comprises Nb + 2n bosons (20| . 
Various features relevant to the shape coexistence phe- 
nomena in the Pb region have been investigated within 
the IBM configuration mixing model: the empirical col- 
lective structures from Po down to Pt isotopes (21rl25| . 
geometry and phases [26U281 ] . and the al geb raic features 
(in terms of the so called intruder spin) |29l43l| . Never- 
theless, one of the main difficulties was that, since too 
many parameters are involved in this prescription, one 
has to specify the form of the model Hamiltonian rather 
a priori and/or simply select the parameters by the fit to 
available data. 

On the other hand, the energy density functional 
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(EDF) framework has been successful in the self- 
consistent mean-field study of bulk nuclear properties 
and collective excitations [32| with various classes of ef- 
fective interactions, e.g., Skyrme [33|,[34[, Gogny [Hj], and 
those used within relativistic mean-field (RMF) models 
[36l I37I ] . The mean-field approximation provides the co- 
existing minima in the potential energy surface, which are 
associated to the different intrinsic geometrical shapes 
[38| . Deformation properties and collective excitations, 
relevant to the shape coexistence in the neutron-deficient 
Pb and Hg isotopes, have already been investigated using 
Skyrme [la 111, Gogny |4lMi^ and the RMF interactions 
[45| as well as within the Nilsson-Strutinsky method j38| . 

Recently a method has been proposed that incorpo- 
rates the merits of the EDF framework into the IBM by 
mapping the potential energy surface (PES) based on a 
given EDF onto the IBM system to generate, energy lev- 
els and transition rates with good symmetry quantum 
numbers [46|. The original idea of Ref. liq has been 
extended to well deformed j47| and triaxial[48[ systems 
and more recently, to the mixing of several multiparticle- 
multihole configurations in Lead isotopes [HI on the basis 
of the Duval-Barrett's technique. 

In this paper we apply the method of Ref. [49( to 
analyze the relevant spectroscopy related to the coexis- 
tence of different intrinsic shapes in the neutron-deficient 
i80-i90jjg i so topes. The predictive power of the model is 
also assessed by comparing its outcome with the available 
data. In addition, the optimal choice of the configuration 
mixing IBM Hamiltonian consistent with the EDF-based 
calculations is identified. Finally, the role of non-axial de- 
grees of freedom in the microscopic EDF PES is stressed. 
We shall use the parametrization DIM of the Gogny-E DF 
[50| . that has been shown (see, for instance (5lT[53l |) to 
have a similar predictive power in the description of nu- 
clear structure phenomena as the more conventional D1S 
[54| parameter set. 

This paper is organized as follows: our theoretical 
framework is briefly summarized in Sec. HU Our results, 
for the microscopic (i.e., EDF) and the mapped PES in 
the considered Hg isotopes, the energy-level systematics, 
the detailed spectroscopy including E2 transitions, the 
spectroscopic quadrupole moment as well as the compar- 
ison with the phenomcnological IBM are given in Sees. 
Mil and [TV] respectively. Finally, Sec. |V] is devoted to 
the concluding remarks. 



II. FRAMEWORK 

We first perform a set of constrained Hartree-Fock- 
Bogoliubov (HFB) calculations using the Gogny-DIM 
[50j | EDF to obtain the corresponding mean-field PES in 
terms of t he g eometrical quadrupole collective variables 
q = (/3, 7) |55j . Note that the PES denotes here the to- 
tal mean-field energy, where neither mass parameter nor 
collective potential is considered explicitly, and that we 
only consider unprotected HFB PES. 



Having the Gogny HFB PES, we subsequently map it 
onto the corresponding IBM PES. In order to treat the 
proton cross-shell excitation, we shall use the proton- 
neutron IBM (IBM-2) because it is more realistic than 
the original version of the IBM (IBM-1), which does not 
distinguish between proton and neutron degrees of free- 
dom. The IBM-2 comprises neutron (proton) s v (s„) 
and d v (d n ) bosons, which reflect the collective pairs of 
valence neutrons (protons) pj| . The number of neu- 
tron (proton) bosons, denoted as N v (N v ), equals half 
the number of the valence neutrons (protons). Since the 
Gogny-DIM PES exhibits two minima in the considered 
Hg nuclei, we take into account up to 2p-2h proton ex- 
citations. The doubly-magic nuclei 164 Pb and 208 Pb are 
regarded as boson vacua (inert cores) for the isotopes 
18 °- 190 Hg. Therefore, N n is fixed, N„ = 1 and 3 for the 
Op-Oh and the 2p-2h configurations, respectively. On the 
other hand, N v varies between 8 and 11. 

The IBM Hamiltonian of the system comprising the 
normal (Op-Ofe) and the 2p-2h configurations is written 

as nnm 

H = V 1 H 1 V 1 +V 3 (H 3 + A)V 3 + H mix , (1) 

where Vi (i = 1,3) stands for the projection operator 
onto the N„ — i configuration space. The operator Hi 
is the Hamiltonian for the configuration with N v = i 
bosons 

Hi = eih d + KiQX".« • QX"- 4 + K ' t L-L+Y^ V P PP> , (2) 

where the first term fid = ^2 p dp ■ d p (p = v or n) stands 
for the (i-boson number operator. The second term in 
Eq. ([2]) represents the quadrupole-quadrupole interaction 
between proton and neutron bosons with Qpf = d p s p + 

spd p +Xp,i[dp x dp]^ being the quadrupole operator. The 
third term is relevant for rotationally deformed systems, 
with L = ylO ^2 p [dp x d p ]^ being the boson angular 
momentum operator. The fourth term on the right-hand 
side (RHS) of Eq. ^ stands for the three-body (cubic) 
boson term between proton and neutron bosons 

Vp P p> = £ < P L JM x 4 x d U {L) ■ [<v x J p x J p] (i) '( 3 ) 

L 

which is identified with the one used in Ref. [48| . For 
each p = v and 7r, there arc five linearly indepen- 
dent combinations in Eq. ©, identified by the value of 
L = 0, 2, 3, 4 and 6. In the present study, we only consider 
the L = 3 term, because its classical limit is proportional 
to cos 2 37 the only term giving rise to a stable triax- 
ial minimum at 7 sa 30° [48|. We have also considered 
the three-body interaction only between neutrons and 
protons since such proton-neutron correlation becomes 
more significant in medium-heavy and heavy nuclei. We 

assume K^J,i = K '^,i = K i > f° r simplicity [48| . 

In Eq. {!]), A represents the energy off-set required to 
excite two protons from the Z = 50 — 82 to the Z = 82 — 
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126 major shells. In the same equation, the term H m i x 
stands for the interaction mixing between the normal and 
the 2p-2h configurations 



(iV,,-,^ , jS, 7) l-Hznixl (iV^-.i, j9, 7)) (i^i'), and reads 



p 3 k4 • 4 + uj d dl • 4)A + h.c. 



(4) 



where the parameters tu s and Ud are the mixing strength 
between the N n = 1 and the N„ = 3 configurations. 

It should be noted that, for each configuration, the 
Hamiltonian in Eq. @ adopts the simplest possible form, 
with a minimal number of parameters consistent with 
the most relevant topology of the EDF PES. Up to the 
third term, the RHS of Eq. ([2]) is the standard form fre- 
quently used in a number of IBM-2 calculations [l6| . In 
the present study, we include the so called L ■ L term since 
we have a relatively large number of bosons, which leads 
to a deformed rotational spectrum. On the other hand, as 
we will see below, the microscopic PES of the considered 
Hg nuclei exhibits a triaxial minimum, which requires the 
inclusion of the cubic term in the IBM Hamiltonian [48[ . 
The physical significance of both the L ■ L and the cubic 
terms has been discussed in detail in Refs. [13] and [48j . 
respectively. 

The geometrical picture of a given IBM Hamilto- 
nian is provided by the coherent-state framework [56j |. 
Such a coherent state represents the boson intrinsic 
wave function specified by the deformation variables 
q = (fa, ^tdTi/jTw)- One can take fa = fa = j3 and 
lu = 7w = 7 a $ only the isoscalar motion is of relevance 
in the following argument. The deformation (3 is assumed 
to be proportional to the one obtained within the HFB 
approximation while 7 is taken the same in the HFB and 
IBM models 0. 

In the configuration mixing model of IBM, one needs 
to consider the direct sum of the coherent state for the 
configuration with N-K.i = i proton bosons [26[, denoted 
here as \§>i(Nnj, fa 7)). The PES is obtained as the lower 
eigenvalue of the 2x2 coherent-state matrix [2(| : 



E(f3,7) 



_(E 11 (fal) E 31 (fa 

EM £33 (A 7)+ A 



(5) 



The diagonal matrix element Ea(fa 7) = 
($ 4 (A^ 4 ,/3, 7)1^1^(^/3, 7)) on the RHS of Eq. © 
is given by 



! + /? 



2 



2 - 2 

4 - 2\j ~(Xu,i + Xn,i)Pi cos 37 + -Xu,iX7t,ift 



K '?N u N Vii (N v + N v4 - 2) Pi sin 2 3 7 , (6) 

(1 + (3f) 3 



where e[ = + and fa = Cp.ifi with 

Gp i being the proportionality coefficient. The non- 
diagonal matrix element is given by Ea> (ft, 7) = 



em = em 



{N^l + 1)^,3 (^j 

1 + hh 
V(i+/3D(i + ^) 



PI 



(7) 



We also assume hereafter w s = ujd — w for simplicity. 

The parameters ej, Kj, x^.ij Xtt,!, k", w and A are 
determined through the same procedure as in Ref. (49j . 
except that, in the present work, the A parameter deter- 
mined from the PES is consistently used for the Hamil- 
tonian [Eq. (JTJ] to be diagonalized. The L ■ L coeffi- 
cient k! is determined independently of the other param- 
eters, that are derived from the PES. It is taken to match 



the Thouless-Valatin moment of inertia [57], calculated 
within the HFB approximation, to the one in the crank- 
ing approximation of the IBM-2 [58[ (i.e., the one that is 
associated to each minimum on the PES [47l |59|). 

Since the oblate HFB minimum occurs always at 
smaller deformation (3 (sa 0.15) than the prolate one 
[ft ~ 0.3) in the considered Hg isotopes, the Hamilto- 
nians for the Op-Oh and 2p-2h configurations are associ- 
ated to the oblate and the prolate minima, respectively. 
On the other hand, the locations of the oblate and the 
prolate minima on the (3 axis remain almost unchanged 
for all the considered nuclei. As a consequence, the scale 
factors Cpj remain almost constant for all nuclei, i.e., 
C/3,i ~ 3 and Cp^ w 5. 

Finally, when the parameters are determined for an in- 
dividual nucleus, the Hamiltonian in Eq. ([T]) is diagonal- 
ized in the enlarged model space in the boson m scheme. 
This gives rise to the energy spectra and wave functions 
for the excited states, that can be used to compute other 
spectral properties. 



III. POTENTIAL ENERGY SURFACES 

The Gogny-DIM and the mapped PES are plotted in 
Fig. [T] as functions of the q = (0, 7) deformations. We 
have restricted the plots to configurations up to 3 MeV 
from the global minimum. Both microscopic and mapped 
PESs give fa values consistent with earlier calculations 
such as the Nilsson-Strutinsky method [38| and the col- 
lective model approach with the Gogny-DlS EDF [4l| . 
where fa ~ —0.15 and 0.25-0.3 for the oblate and pro- 
late configurations, respectively. 

The Gogny-DIM PESs for !80,i82 Hg display notable 
prolate minima at j3 ~ 0.3. On the other hand, the 
oblate minimum at (3 rj 0.15 becomes gradually lower 
in energy when approaching 184 Hg, whose PES becomes 
softer along the 7 direction. The prolate-oblate energy 
difference in the HFB PES reaches a minimum in 186 Hg 
which, at least within the Gogny-DIM HFB framework, 
represents the most dramatic example of shape coexis- 
tence in the considered Hg chain. Note that the lower 
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FIG. 1: (Color online) Microscopic ("DIM") and mapped ("Mapped") potential energy surfaces for the isotopes 180 190 Hg 
in the (/3, 7)-plane are plotted up to 3 MeV from the absolute minimum. The microscopic results are obtained with the 
Gogny-DIM EDF. 




FIG. 2: (Color online) Derived IBM parameters (a) ti, (b) Ki, (c) (d) Xir,i> ( e ) K i > (f) K 'ii (g) A an d (h) w f° r the considered 
180-190 jjg nuc i e j ag f unc tions of mass number A. Note that, in panel (f), the parameter n'i is plotted in keV unit. Figure legends 
in panels (a-f) are indicated in panel (a). 



minimum becomes a bit off from the prolate axis with 7 ss 5°. In 188 Hg, one can clearly sec that the oblate 
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minimum becomes energetically favoured over the one 
around 7 sa 10° on the prolate side. In 190 Hg only the 
oblate minimum survives. The mapped IBM PES nicely 
follows the topology of the Gogny-DIM PES (i.e., loca- 
tions of the minima as well as the curvatures around a 
given minimum in both the /3 and 7 directions) for an 
individual nucleus. The overall systematic trend with 
neutron number is also well reproduced. 

The derived IBM parameters [Eq. (JTJ] are depicted in 
Fig. [5] as functions of the mass number A. In agreement 
with its empirical boson number dependence |l6[ , as well 
as with our previous findings (4(|, the single d— boson 
energy e shown in panel a) exhibits a parabolic behavior 
with respect to the mid-shell. This also agrees, with the 
empirical evolution of the 2 + excited state expected in a 
given isotope/isotone sequence. As can be observed, the 
ei and £3 parameters, shown in the figure, roughly follow 
this empirical trend. Contrary to the fitting calculations 
[20I |2lj . the parameter €3 is always larger than e\. The 
parameters fti.3, shown in panel b), do not change too 
much. Nevertheless, they are several times larger than 
the phenomenological ones ( k,\ ss —0.17 ~ —0.14 MeV 
and « 3 w -0.14 - -0.11 MeV) [22|. This is a direct con- 
sequence of the topology of the Gogny-DIM PES, whose 
energy minimum is too deep leading to a large curvature 
around the minimum on the j3 axis. 

Generally, whether a given nucleus is prolate or oblate 
is specified by the sign of the sum x„ + x-k- As seen 
from Figs. [^c) andj^d), the sum has positive (negative) 
sign for the oblate (prolate) configurations, being consis- 
tent with the microscopic PES. In a previous study on 
Hg isotopes (e.g., [Hjj), the 2p-2h configuration has been 
considered in the rotational SU(3) limit of the IBM [l6[ 
by taking x = V</2 ~ 1.3. The present result seems to 
contradict this assumption since both Xu an d Xn values 
for the 2p-2h configuration are away from the SU(3) limit 
of V7/2 reflecting the 7 softness of the Gogny-DIM PES. 
Thus, not only the x parameters but also the cubic term 
should be important for the description of the triaxiality 
in the considered Hg isotopes. 

From Fig. [21(e), one sees that the derived k'{ value for 
both 186 Hg and 188 Hg is particularly large in agreement 
with the Gogny-DIM PES [FigH] of the two nuclei dis- 
playing the most notable 7 softness in the considered 
isotopic chain. The L ■ L coefficient k', shown in panel f), 
appears to be almost constant, but its sign is of partic- 
ular relevance here. The positive (negative) sign for the 
normal (2p-2h) configuration implies that the inclusion of 
the L ■ L term reduces (enlarges) the moment of inertia of 
the unperturbed collective band of the configuration. As 
we will show later fSec. IIV|) . the L ■ L may have a crucial 
impact on the level spacing for each band. 

The energy offset A in Fig. [21(g) changes with neutron 
number symmetrically with respect to 184 Hg. The energy 
needed to excite two protons across the Z = 82 closed 
shell becomes maximal for this mid-shell nucleus because 
the intruder 2p-2h configuration gains maximal energy 
through deformation. Finally, as can be observed from 
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FIG. 3: (Color online) Level-energy systematics for 180_190 Hg 
isotopes with mass number. Theoretical level energies re- 
sulting from the mapped IBM-2 Hamiltonian based on the 
Gogny-DIM EDF (a) are compared with the experimental 
(b) energies taken from the ENSDF database [gfj] and from 
Ref. [T2| ■ To guide the eye, dashed and solid lines connect the 
yrast and the non-yrast states, respectively. 



panel h) , the mixing strength uj remains rather constant 



IV. SPECTRAL PROPERTIES 

Let us now discuss how the topology, as well as the 
evolution, of the microscopic PES and the derived pa- 
rameters correspond to the spectroscopic properties of 
the considered Hg isotopes. 



A. Level energy systematics 

We show in Fig. [3] the low-lying excited states of the 
isotopes 18 °- 190 Hg with spins J < 6 + . The theoretical 
energy levels in Fig. [3ja) , obtained through the diago- 
nalization of the mapped IBM-2 Hamiltonian, are com- 
pared with the corresponding experimental data, shown 
in panel b) , taken from the Evaluated Nuclear Structure 
Data File (ENSDF) (6(J and from Ref. [3 for 180 Hg. 

We observe, a reasonable qualitative agreement be- 
tween the calculated and the experimental levels. The 
lowest two + states in the neutron-deficient Hg isotopes 
originate cither from the oblate or the prolate configura- 
tions. Starting from 190 Hg down to the mid-shell nucleus 
184 Hg, the calculated yrast states come down following 
the trend of the experimental levels. A rather good agree- 
ment with the experimental data has been obtained for 
Hg nuclei with A > 184, including the ordering of the 
energy levels and the location of the 0^" state. Experi- 
mentally, the 0^ excited state in 18 °. 182 Hg (the bandhcad 
of the intruder prolate band) is still lying close to, and 
almost degenerate with, the 2^~ excited state around the 
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FIG. 4: (Color online) Fraction of the intruder 2p-2h config- 
uration in the wave functions of the 0f 2 i the 2f 2 an d the 
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excited states of the considered 



Hg isotopes (in 



%). The bars for each nucleus (from left to right) correspond 
to Of , Of, 2J~, 2f , 4f and 4f excited states, respectively, as 
indicated in the figure. 
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FIG. 5: (Color online) Level scheme for the 188 Hg nucleus. 
The level energies and the B(E2) values (numbers written 
along arrows) obtained from the mapped IBM-2 Hamilto- 
nian based on the Gogny-DIM EDF are compared with the 
corresponding experimental data The oblate and the 

prolate bands are indicated by "Obi." and "Prol.", respec- 
tively. To help identify the corresponding values in the plot, 



B(E2;4+ -> 2+) = 84 W.u. and B(E2;0^ 
Weisskopf units) in the right panel. 
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middle of the major shell. On the other hand, the theo- 
retical non-yrast levels increase their excitation energies 
for A < 182. In particular, the Of excitation energy is 
too large in comparison with the corresponding experi- 
mental value. This can be due to the fact that the mean- 
field oblate-prolate energy difference predicted with the 
Gogny-DIM EDF could be too large for these two nu- 
clei. In general, the HFB oblate-prolate energy differ- 
ence should roughly correspond to the excitation energy 
of the Of before mixing. Once the two configurations are 
mixed, the + excitation energy becomes more or less 
larger than the energy difference between the two mean- 
field minima. This is due to the level repulsion effect 
that is seen in 186 > 188 Hg. However, the Of excitation en- 
ergy for 18 °d82,i84-r_jg seems t ]-, e a l mos t equal to, and 

even smaller than, the oblate-prolate energy difference 
in the corresponding PES. The reason is probably that, 
mainly due to the triaxiality and to the shift of the d bo- 
son energy by 6k' necessary to keep the topology of the 
PES, the unperturbed Of excitation energy is sufficiently 
smaller than the oblate-prolate energy difference in the 
HFB PES. As a result, the energy of the Of level after 
the configuration mixing, may be even smaller than the 
energy difference between the two mean-field minima. 

In particular, to assess the effect of the cubic term on 
the spectrum, we performed a test calculation, in which 
the cubic term was not included in both configurations. 
Without the cubic term, the derived mixing strength u> 
turned out to be several tens per cent larger than the 
one in Fig. [5{h) in order to reproduce the topology of 
the EDF PES along the 7 direction. This resulted in 
stronger level repulsion between the two + eigenstates, 
and thus the Of excitation energy was overestimated by 
several hundred kcV in comparison with Fig. E^a). 



B. Detailed level scheme 

To identify band structures, we now analyze the de- 
tailed level scheme of individual nuclei, including both 
the inter- and the intra-band E2 transition rates. We 
will consider the isotopes 182 > 18 4,iS6,i88jjg near the neu- 
tron mid-shell N rj 104. These isotopes exhibit, the most 
spectacular competition between different shapes among 
the considered Hg nuclei. To facilitate the comparison 
with the experimental data, the excited states shown be- 
low are classified into either oblate or prolate band ac- 
cording to the prolate-oblate predominance of the wave 
function of the state and/or the E2 transition sequence. 

To this end, we have first computed the probability 
of the basis states in the wave function of the state of 
interest. In Fig.2]we have plotted, the fraction of the 2p- 
2h component in the wave functions of the 0f 2 ; 2f 2 an d 
the 4f 2 states. One sees that the yrast states in 18 °. 182 Hg 
are composed of the intruder prolate configuration. The 
two configurations are mixed for each of the low-lying 
states of 184 > 18 6jjg On the other hand, for the nucleus 
188 Hg the tendency is the opposite to the one in 182 Hg. 
Only a small fraction of the intruder component plays 
a role in the low-lying states of 190 Hg. It appears that 
the excited + state of the nucleus 190 Hg is originated 
from a single oblate configuration, consistently with the 
empirical observation |4lj . 

The -B(E2; J — > J') transition rate reads 

B(E2; J -> J') = _L_|(J'||f( E2 )||J)| 2 , (8) 

where | J) (| J')) represents the wave function of the ini- 
tial (final) state with spin J (J'). The E2 operator T < E2 ' 
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is given as f ^ = E P , 4 =i, 3 e p ^V t Q x P ^P u where Q* pi is 
identified with the quadrupole operator in Eq. ([2]). We 
considered the same \p,i value as the one used in the diag- 
onalization of the Hamiltonian in Eq. ([TJ . This is equiva- 
lent to the so-called consistent-Q formalism in the IBM-1 
foil ]. The boson effective charge e Pi j is assumed to be the 
same for protons and neutrons, i.e., e v ^ = e w> i = The 
ei value is fixed to e\ = 0.166 eb, taken from Ref. f22| . 
while the e% value is determined for each individual nu- 
cleus according to the empirical relation e^/ei ~ I/C3/K1I 

We show in Fig. [5] the detailed level scheme of the 
nucleus 188 Hg. The spectra and the £?(E2) transition 
strengths, computed from the mapped IBM-2 Hamilto- 
nian based on the Gogny-DIM EDF, are compared with 
the experimental data in the same figure. A nice agree- 
ment between theory and experiment is obtained regard- 
ing the band structure, the energy of the 0^" state as 
well as the level ordering. A pronounced mixing be- 
tween the two configurations can be confirmed from the 
_B(E2) values for the transitions 6^~ 2 — > 4+ 2 . They re- 
flect the sizable amount of the mixing seen in the wave 
functions of the two 6 + excited states [the wave func- 
tion of the 6^ (63 ) state consists of 69.1 (38.2)% of the 
2p-2h component]. Mainly due to the uncertainty in the 
boson effective charges, the i?(E2) values for the transi- 
tions 2\ — > 2\ and 2\ — > Of are rather larger than the 
corresponding available data. On the other hand, the 
branching ratio for these transitions compare well with 
the experiment. 

Consistent with the experimental evidence, the + 
ground-state in 186 Hg (see, Fig. EJ) is predicted to be of 
oblate nature. The normal and the intruder configura- 
tions are substantially mixed in the wave functions of the 



lowest two + states. According to Fig. |H the 0^ (0^ ) 
state comprises 36.7 (62.3)% of the intruder 2p-2h con- 
figuration. Actually, the prolate configuration is quite 
notable in the mean-field PES shown in Fig. [TJ In addi- 
tion, our calculations seem to reproduce the experimen- 
tal level scheme including the £>(E2) pattern: the oblate 
ground-state band is rather stretched, while the intruder 
prolate band built on the excited state is compressed. 

Such systematics in the level spacing for the lower band 
in each configuration, which is likely to be the K 71 = + 
band, can only be reproduced with the inclusion of the 
L ■ L term. Actually, we have also performed a test calcu- 
lation for all nuclei considered, that did not include the 
L-L term, and especially the intruder prolate band turned 
out to be, as expected, too stretched inconsistently with 
the data. This finding represents, of course, only a prac- 
tical side of the L ■ L term. More importantly, this result 
implies that to describe the shape coexistence and hence 
the configuration mixing in the considered Hg isotopes 
within the present model, not only the static HFB PES 
but a specific dynamics should be also taken into account 
through the L-L term with nonzero angular frequency. In 
addition, a pronounced mixing between the two configu- 
rations is observed from the branching ratios of the E2 
transition, such as for the 4^~ 2 — > 2^ 2 transitions, which 
are indeed comparable with each other. 

We also address the study of collective bands, other 
than those built on the two + states, for 186 Hg as there 
are sufficient experimental data to compare with. From 
Fig. [5] one sees that the present calculation also repro- 
duces the excitation energies with odd spin quite nicely. 
The odd-spin states, 3f, 5+ and 7± (3^, 5^ and 7 2 ), 
are predicted to be the members of the prolate (oblate) 
band. Due to the strong E2 transition sequence, it is 
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FIG. 7: (Color online) Same as Fig. [5] but for 184 Hg nucleus and B(E2; 3+ -> 2+) = 153 W.u. 



also likely that a set of the states 2% , 3f, 4^", 5+ , 63" and 
7 1 (24 , 3^, 44 , 5^, 64 and 7^), forms a quasi-7, i.e., 
K"* = 2 + , band for the prolate (oblate) configuration. 
These two quasi-7 bands seem to be close in energy, with 
2 + bandheads being within 400 keV. One can also observe 
in both quasi-7 bands quite strong E2 transitions from 
the 3 + level to the corresponding 2 + bandhcad, which 
is in the same order of magnitude as the 2 + — s- + E2 
transition in each /(" = + band, and those between 
the members of the quasi-7 band. Note that this pre- 
diction (i.e., the existence of the two quasi-7 bands) is 
consistent with the empirical assignment of these levels 
[4l| . including the collective model description based on 
Gogny-DlS EDF. In addition, one should notice that the 
quasi-7 band in the prolate configuration looks similar 
to the one predicted within the rigid-triaxial rotor model 
of Davydov and Filippov, that is, the doubles (2+ ,3+), 
(4+,5+), (6+,7+), ... etc (62|. From the comparison with 
the IBM phenomenology for 186 Hg [22 |. one sees that 
the present result exhibits a similar level of agreement 
with the experimental one regarding the energy spectra 
of the oblate ground-state band with J < 6 + . Our result 
also better reproduces the 3+ state while, as addressed 
in Sec. IIV A[ the prolate 0^ band-head energy is a bit 
more overestimated. 

From Fig. [71 one observes that the quality of agree- 
ment with the experiment in 184 Hg is similar to the one 
obtained for 186 Hg. From the strong E2 transitions to 
the corresponding + bandheads, we assume that the 2+ 
and the 2j states belong to the oblate and the prolate 
bands, respectively. As seen from Fig. [H although the 
wave function of the 2+ (2 J) state is predominantly of 
the 2p-2h (Op-Oh) nature, the mixing between the two 
configurations is quite strong: 65.4 (37.5) % of the 2f 
(2 J) wave function is comprised of the intruder 2p-2h 
configuration. Again, the strong mixing between the two 
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configurations is confirmed by the branching ratios of the 
4^2 to the 2+2 excited states. Similar to 186 Hg, we ob- 
serve a significant role played by the L ■ L term in the 
description of the level spacing. As seen in Fig. [71 the 
present study predicts the existence of the two quasi-7 
bands for 184 Hg as well. It should be of interest to iden- 
tify such band structure experimentally. In the present 
framework, such systematics of the 7 band seems to be 
seen most vividly in 184 > 186 jjg since the oblate and/or the 
prolate mean-field minimum is the softest and the closest 
in energy among the considered Hg isotopes. Regarding 
the nucleus 188 Hg, however, any sequence of the quasi-7 
band structure has not been obtained in our calculations 
(not shown in Fig. [S]). This agrees well with available 
data [6(|. On the other hand, the predicted 3+ en ergy 
(1.437 MeV) agrees well with the data (1.455 MeV) (6Qj- 
In Fig. [51 we show the results corresponding to the nu- 
cleus 182 Hg. Rather contrary to the experimental assign- 
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moment for a state with spin J reads 

& S) ^ = \IW(-J o j)( J W fiE2) \\ J )- 



(9) 



The calculated quadrupole moments for the 2f and the 
2% states in the 18 °- 190 Hg nuclei, are shown in Fig. |H1 
For the lighter isotopes !80-i84 H& Q(«)(2+) [Q^(2+)j 
is negative (positive), reflecting the prolate (oblate) na- 
ture of the state. On the other hand, Q^(2^) becomes 
positive for A > 186. For both i84,i86 Hg) the predicted 
Q( s ) values are smaller in magnitude than in other nu- 
clei, as there is significant triaxiality in these nuclei. The 
overall systematic trend in Q^(2^ 2 ) seems to correlate 
well with the evolution of mean-field minima shown in 
Fig. [T] and with the structure of the corresponding wave 
functions in Fig. 2J 

To compare directly with the experimental data, one 
can also deduce the intrinsic deformation parameter 
foiJ) from Q( s >(J), through the relation 



ment, the present calculation predicts the ground-state 
band to have an intruder prolate nature, as 75.4% of the 
0^ state is dominated by the 2p-2h configuration (see, 
Fig. [4j. A clear collective pattern is seen from the cal- 
culated E2 transition sequence for both the prolate and 
the oblate bands, although the B(E2;2^ — > 0+) value 
is rather overestimated. Another major deviation from 
the experimental data occurs in the intra-band transition 
from the 4^ state in the oblate band to the 2+ state in 
the prolate band. The corresponding experimental value, 
B(E2;4^ 2f) = 176(19) W.u, is quite large in com- 
parison to the 2^ — > E2 transition strength reflect- 
ing the very strong mixing between the different config- 
urations. In our calculations the intra-band transition 



S(E2;4+ 2+) = 2 W.u., while B{E2 



2+ 
z i 



0J) = 116 



W.u., although there exists a certain amount of mix- 
ing between the two configurations in the corresponding 
eigenstates (see, Fig. 2]). A possible reason for this devi- 
ation could be that, in comparison to the experimental 
data, the theoretical 0^" excitation energy is too large, 
leading to a weak E2 transition. This problem arises 
mainly because the energy difference between the oblate 
and the prolate minima is too large in the HFB PES 
(see, Fig. [T|). A similar argument holds for 180 Hg, whose 
wave functions for the low-lying states look very similar 
in the composition of the two configurations to the ones 
for 182 Hg in Fig. |4l exception made of the 4+ state. 



C. Spectroscopic quadrupole moment 

In order to better understand the deformation prop- 
erties of the considered Hg isotopes, we have also ana- 
lyzed the spectroscopic quadrupole moments for the low- 
est two 2 + excited states belonging either to the first 
oblate or prolate band. The spectroscopic quadrupole 



&(J) 



5 4ttQW(J) 
Tott 3ZR 2 



(10) 



where R = 1.2A 1 / 3 fm 2 and Q^(J) represents the 
trinsic quadrupole moment given by 



m- 



QW(J) 



2J + 3 



>(«) 



(J) 



(11) 



with the assumption K = 0. In the present work, the 

calculated /3f w (2 + ) value is generally consistent with the 

one taken from the mean-field PES in Fig.[T](/32 ~ 0.15), 

whereas the 02° one is considerably smaller. In 
fact, for i80,i82,i88,i90 Hgj we obta . n _ Q 1? ^ pM{2+) < 

-0.15 and 0.15 < ^ rol {2+) < 0.16, while for 184 > 186 Hg 
both pl ro1 and |/3 2 oW | are much smaller than 0.1. These 
small ,02 (2 + ) values seem to be rather reasonable because 
the effects of the configuration mixing and the non-axial 
degrees of freedom (characterized by the rather small 
magnitudes of and Xw)j which we have emphasized 
so far, are taken into account through the diagonaliza- 
tion. 



D. Comparison with the phenomenological fit 

Finally, we assess the validity of the set of IBM param- 
eters derived here by comparing it to the one in Ref. pl| 
determined from a fit to the experimental spectra. Both 
sets of parameters lead to reasonable agreement with the 
data, but they differ in their description of the ground- 
state shape. 

To illustrate this, we show in Fig. [TU] the PES for 
the nucleus 184 Hg using the set of parameters used in 
Ref. [21 L The IBM Hamiltonian for both configurations 
in Ref. 2M was identical to the one on the RHS of Eq. @ 
up to the second term, but included the neutron-boson 
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FIG. 10: (Color online) The potential energy surface of the 
184 Hg nucleus, based on the parameters determined from a 
phenomenological fit [5l|] to the experimental spectra. The 
PES is drawn in terms of the boson quadrupole deformation 
variables q = 



interaction term V uv = \ J2l=o 2 4 ^1(^1 x ' {d v x 

d v ) {L) . Note that the PES in Figgis plotted using the 
same energy scale as in Fig. Q] but in terms of the boson 
variables q = 7), as they are not connected to the 
geometrical deformation variables. 

One sees in Fig.[10]the spherical minimum, which leads 
to the vibrational U(5) [l6| spectra, as well as the sec- 
ond minimum on the prolate axis, corresponding to the 
intruder deformed band equivalent to the SU(3) limit of 
the IBM pjj]. Therefore, the corresponding energy spec- 
trum can be considered as the coupling of the two limits. 
However, this picture seems to be rather simplistic. If 
one compares the mapped PES in Fig. [1] and the one 
in Fig. [TU] a difference is observed in the role of triax- 
ility and the pronounced oblate minimum at j3 ~ 0.15 
of the mapped PES. Indeed, the mapped PES is gener- 
ally rather sharp in both the oblate and prolate axes and 
exhibits a richer structure along the 7 direction. On the 
other hand, the phenomenological PES appears to be flat 
in both (3 and 7 directions. 

A possible reason for the difference would be the lim- 
ited number and/or species of bosons in the IBM. Re- 
garding the former, an inert core is assumed in the bo- 
son intrinsic (coherent) state and so the boson number is 
already a good quantum number. On the other hand, 
within the self-consistent EDF calculation carried out 
here the model space consists of all nucleons but the par- 
ticle number is still not a good quantum number. In 
the mapped PES, such effect may be renormalized into 
the parameters in the configuration-mixing IBM Hamil- 
tonian through the calibration with the EDF PES and 
so does not matter too much in the present framework. 
In this respect, it should be also interesting to see how 
qualitatively /quantitatively could a symmetry- projected 
calculation change the topology of the mean-field PES 
and the resulting spectrum in the considered Hg nuclei. 
V. SUMMARY 



IBM- 2 with configuration mixing based on the Gogny- 
D1M EDF. The parameters of the IBM Hamiltonian were 
determined through the mapping of the Gogny-DIM PES 
onto the one provided by an IBM-2 Hamiltonian in the 
boson condensate. 

The importance of the 7 degree of freedom has been 
emphasized in the description of the ground-state shape 
of the considered Hg nuclei with the Gogny-DIM EDF. 
Spectroscopic quantities that help to identify the oblate- 
prolate shape coexistence, resulting from the mapped 
Hamiltonian, have been compared with the experimen- 
tal data. The mapped Hamiltonian provides us with 
a fairly good description of the observed energy lev- 
els without any adjustment, including the predominance 
of the oblate or the prolate configuration in each col- 
lective band. Particularly in the near mid-shell nu- 
clei 184 . 186 > 188 Hg, both the calculated oblate and prolate 
bands indicate a notable collective pattern in the in-band 
and the intra-band E2 transition rates, the latter being 
a strong evidence for the mixing between the two config- 
urations. 

It has been shown that the observed collective band 
structure for both the oblate and the prolate configura- 
tions is reproduced only with the proposed Hamiltonian 
including the cubic term and the L ■ L term. The cubic 
term has been introduced following the indications of the 
EDF PES, which reveals the importance of the triaxial 
degrees of freedom. It is also necessary not only to use 
the static HFB PES information, but also to consider the 
response of the system to nonzero angular frequency as 
required by the L ■ L term. In this respect, the present 
work suggests a possible optimal choice of the configura- 
tion mixing IBM Hamiltonian as well as the necessity to 
take into account the non-axial degrees of freedom within 
an EDF framework for the neutron-deficient Hg isotopes. 

Nevertheless, the disagreement between our results and 
the data in 180 ' 182 Hg should not be overlooked. This 
problem has been mainly attributed to the paramctriza- 
tion and/or the form of the considered EDF, which pre- 
dicts that the two HFB minima differ too much in energy. 
An extension of the present EDF framework would help 
to clarify this problem. 
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